home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / taylor.f < prev    next >
Text File  |  2002-06-11  |  598b  |  22 lines
  1. C
  2. C PAGE 307: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: TAYLOR.FOR
  5. C
  6. C TAYLOR SERIES METHOD (ORDER 4) FOR SOLVING AN ORDINARY DIFFERENTIAL EQUATION
  7. C
  8.       DATA  T,X/2*0.0/, H/7.8125E-3/  
  9.       PRINT 3,T,X 
  10.       DO 2 K=1,128
  11.         X1 = 1.0 + X*X + T**3 
  12.         X2 = 2.0*X*X1 + 3.0*T*T       
  13.         X3 = 2.0*X*X2 + 2.0*X1*X1 + 6.0*T       
  14.         X4 = 2.0*X*X3 + 6.0*X1*X2 + 6.0 
  15.         X = X + H*(X1 + H*(X2/2.0 + H*(X3/6.0 + H*X4/24.0)))
  16.         T = REAL(K)*H       
  17.       PRINT 3,T,X 
  18.    2  CONTINUE
  19.    3  FORMAT(5X,F10.5,5X,E22.14)      
  20.       STOP
  21.       END 
  22.